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We derive the effective low-energy theory for interacting 
electrons in metallic single- wall carbon nanotubes taking into 
account acoustic phonon exchange within a continuum elastic 
description. In many cases, the nanotube can be described as 
a standard Luttinger liquid with possibly attractive interac- 
tions. We predict surprisingly strong attractive interactions 
for thin nanotubes. Once the tube radius reaches a critical 
value _Ro ~ 3.6 ± 1.4 A, the Wentzel-Bardeen singularity is 
approached, accompanied by strong superconducting fluctua- 
tions. The surprisingly large -Ro indicates that this singularity 
could be reached experimentally. We also discuss the condi- 
tions for a Peierls transition due to acoustic phonons. 

PACS number(s): 73.63.Fg, 72.10.-d, 74.25.Kc 



I. INTRODUCTION 

Superconductivity in carbon nanotubes has lately gen- 
erated a lot of excitement. Experimental observations 
include an unexpectedly strong proximity effect in single- 
wall nanotubes (SWNTs) [1,2], fluctuation supercon- 
ductivity in ultrathin SWNTs [3], and intrinsic BCS- 
type behavior in "ropes" containing many SWNTs [4]. 
These experiments have in turn led to a number of 
theoretical papers [5-8]. In particular, the authors of 
Refs. [5,6] argue that phonon exchange holds respon- 
sible for some attractive electron-electron interaction. 
The value or even order of magnitude of the electron- 
phonon coupling was however left open or simply ex- 
tracted by fitting experimental curves. In this work, we 
wish to fill this gap and provide a quantitative analy- 
sis of phonon-mediated retarded electron-electron inter- 
actions in a metallic SWNT. (The theory also holds qual- 
itatively for a heavily doped semiconducting SWNT.) 
Below we include the usual (repulsive) Coulomb inter- 
actions in the effective low-energy theory valid at ener- 
gies \E\ < vp/R, with Fermi velocity vp and SWNT 
radius R. We assume that the SWNT is doped away 
from the charge neutrality point, Ep ^ 0, i.e., an in- 
commensurate situation, which is normally encountered 
in practice. Then electron-phonon and electron-electron 
Umklapp processes can be neglected. 

Besides mediating possibly attractive interactions, the 
effect of phonons on the electronic system could in princi- 
ple lead to two instabilities, namely the Peierls transition 
[9] and the Wentzel-Bardeen (WB) singularity [10]. All 



these issues can be investigated by integrating out the 
phonons and studying the electron system alone, since 
the relevant backscattering electron-phonon interactions 
are weak [11]. This is the approach taken in our paper. 
Our treatment of phonons expands on recent work by 
Suzuura and Ando [12], where the field theory for acous- 
tic phonons and their coupling to electrons in SWNTs 
was derived. Optical phonons are ignored in this theory 
as well as in ours. This should create no major prob- 
lem for the issues at stake here, since optical phonons 
do not produce sizeable attractive interactions nor a WB 
singularity. Their relevance for the Peierls transition is 
discussed briefly below. The theory of Ref. [12] is ana- 
lytically tractable and, moreover, appears to be in excel- 
lent agreement with local-density-functional calculations 
even for ultrathin (i? w 2 A) nanotubes [13]. For in- 
stance, we find that the frequency of the breathing mode 
computed in Ref. [13] coincides to within 4% with the 
value extracted from the continuum theory of Ref. [12]. 
Furthermore, the predicted value for the twist velocity 
vt below is consistent with Refs. [14,15]. This is for- 
tunate since previous estimates for the electron-phonon 
coupling [16,17] were only qualitative in nature and did 
not even agree on orders of magnitude. For related the- 
oretical work on acoustic phonons and electron-phonon 
interactions in SWNTs, see Ref. [18]. 

As we show in detail below, the phonon-mediated ef- 
fective attraction among electrons is very strong for suffi- 
ciently thin nanotubes. The dominant phonon exchange 
processes involve two phonon modes, namely the stretch- 
ing and the breathing mode. The stretching mode has a 
linear dispersion characterized by a velocity vs, while the 
breathing mode has a finite frequency lob at long wave- 
lengths. Following general arguments [11,19], for a large 
part of the relevant parameter space, the influence of the 
stretching mode is negligible. Then the system can be 
described as a Luttinger liquid (LL) with standard inter- 
action parameter Kc+ in the charged channel [20] , where 
Kc+ — 1 for noninteracting electrons. Phonon exchange 
involving the breathing mode implies 
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(1.1) 



with Rb 2.4 ± 0.9 A. Here K^^ < 1 characterizes 
the repulsive Coulomb interactions, where K^_f_ ~ 0.2 
for unscreened interactions [20]. Notice that without ex- 
ternal screening, phonon exchange is not relevant, and 
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Kc-\- ~ J^c+- However, once Coulomb interactions are ex- 
ternally screened off, e.g. by nearby gate electrodes, we 
get K^_^ Ri 1 and Eq. (1.1) predicts surprisingly strong 
attractive interactions {Kc+ > 1). For a (10,10) arm- 
chair tube with R = 6.79 A, the LL parameter is then 
Kc+ ~ 1.3, i.e. significantly larger than the noninteract- 
ing value. These attractive interactions get stronger for 
thinner nanotubes. 

As we discuss below, the neglect of the retarded in- 
teraction due to the stretching mode, leading to a (four- 
channel) LL with interaction parameter Kc+ given by 
Eq. (1.1) and interaction parameters K\ = 1 in the 
three neutral channels (sec below), is valid provided R 
is sufficiently large compared to the "critical radius" 
Rc = {K^^)'^Ro, where Rq is slightly larger than Rb, 
Rq = 3.6 ± I.4A. Only in the immediate vicinity of this 
radius, the retarded interaction mediated by the stretch- 
ing phonon mode has to be kept explicitly. This criti- 
cal radius defines the location of the so-called Wcntzcl- 
Bardeen (WB) singularity [10]. For R < Rc, our model 
breaks down, and a transition to a phase-separated state 
is expected. Previous work by Loss and Martin [21] stud- 
ied the WB singularity in detail for one-dimensional (ID) 
systems, but a concrete physical realization seemed out 
of reach due to the small ratio between sound and Fermi 
velocity in all systems considered so far. We argue be- 
low that SWNTs with screened-off interactions may of- 
fer the possibility to reach this elusive singularity, which 
has never been observed experimentally. This is remark- 
able since the ratio between sound and Fermi velocity 
is still small. However, in a SWNT the stretching and 
the breathing phonon modes collaborate in driving the 
system towards the WB singularity in a very efficient 
manner. 

Moreover, the closeness to the WB singularity holds re- 
sponsible for rather strong superconducting fluctuations. 
From our theory, estimates for the attractive interactions 
generated via phonon exchange are obtained. To de- 
scribe most of the issues raised by the experiments of 
Rcfs. [3,4] in a quantitative way, it will however be im- 
portant to also analyze the effect of various inter-tube 
couplings in SWNT ropes or crystals in more detail. Em- 
ploying a mean-field approximation for Josephson cou- 
plings, the standard ID Ginzburg-Landau description 
can be obtained [22]. In this paper, we focus on the 
phonon-mediated electron-electron interaction within a 
given nanotube, and postpone a full discussion of exper- 
imentally relevant quantities to a future publication [23]. 

The present paper also sheds light on the recent dis- 
cussion about the possibility of a Peierls transition in 
SWNTs [15,24-26]. Notably, even in the incommensu- 
rate situation studied here, in principle a Peierls transi- 
tion can occur [27], although the resulting pscudogap is 
expected to be wiped out by quantum fluctuations [11]. 
Within a mean-field approximation, the Peierls transi- 



tion due to optical phonons has been studied in Refs. 
[24-26], resulting in a Kekule-type modulation. Since 
we only treat acoustic phonons, we have nothing to say 
about this issue. However, recent theoretical work includ- 
ing acoustic phonons has pointed out the possibility of a 
sizeable Peierls gap in thin armchair SWNTs [15]. This 
Peierls modiilation was found to correspond to a twist 
distortion of the SWNT. Below we study arbitrary metal- 
lic SWNTs, and indeed confirm that the relevant Peierls 
distortion for an armchair SWNT would correspond to a 
tube twisting. In addition, we find the proper mode for 
any metallic SWNT. Our findings follow from analyzing 
backscattering phonon exchange, which operates at the 
doping-dependent wavevector 2qp, where qp = \Ef\/vf- 
Based on our essentially exact results, we do not expect 
that the Peierls transition predicted in Rcf. [15] is ob- 
servable, at least not for an individual SWNT away from 
the neutrality point. According to our theory, a Peierls 
gap in the ID sense [11] could only be possible for very 
thin SWNTs and under rather stringent conditions. 

The outline of the paper is as follows. In Sec. II 
we briefly summarize the elastic contimiTim theory for 
phonons and the electron-phonon interaction. The re- 
sulting effective electron-electron interaction generated 
by phonon exchange is derived in Sec. III. The domi- 
nant forward scattering process, the Luttinger liquid with 
attractive interactions, and the Wentzel-Bardeen singu- 
larity are discussed in Sec. IV. The remaining phonon 
exchange terms and the Peierls transition is studied in 
Sec. V. Some details concerning Sec. V have been trans- 
ferred to an appendix. Finally, we discuss the relevance 
of inter-tube phonon exchange and conclude in Sec. VI. 

II. ACOUSTIC PHONONS AND 
ELECTRON-PHONON COUPLING 

Here we shall always assume a suspended (free- 
standing) tube. With regard to the important phonon 
modes coupling to the low-energy electrons, this should 
be an accurate assumption both for SWNTs embedded 
in a zeolite matrix [3] or in a rope [4]. However, modifi- 
cations may arise for individual SWNTs on a substrate, 
where phonon modes could be pinned. Let us then briefly 
summarize the main results of Ref. [12], where acoustic 
phonons in SWNTs are described by an elastic contin- 
uum theory. The Euclidean action is (we put h = 1) 

S = J drcPr {(M/2)u2 + U[u] + Uc[u]} , (2.1) 

where u^{r) = {ux,Uy,Uz) is the displacement field and 
r = {x,y). The displacement field in this approxima- 
tion does not discriminate among the graphite sublat- 
tices, and hence will not describe optical phonons. The 
coordinate system is chosen as follows. The y-axis points 
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along the tube, < a; < 2'kR around the circumference, 
and z denotes the direction perpendicular to the tube 
surface. The elastic potential energy density in Eq. (2.1) 
is given by 

UM = \ {B{UxX + Uyyf + ll[{U:oX - Uyyf + 4U^y] ) , 

where the strain tensor is defined by the relations Uyy = 
dyUy, 2uxy = dyUx + dxUy, Aud , due to the cylindrical ge- 
ometry, Uxx = dxUx + Uz/R- Furthermore, the curvature- 
induced potential energy density is 

In the above expressions, a = 2.46 A is the lattice 
constant, and the carbon mass per unit area is M = 
3.80 X 10"'^ kg/m^. The bulk modulus B and the shear 
modulus fj. are accurately known for bulk graphite, 



B/M = 2.90 X 10^mVsec^ 
^/M = u| = 1.51 X lO^mVsec^ 



(2.2) 



Finally, the normal-force constant is E/M = 6.19 x 
lO^m^/sec^. Since S <C -B,/i and the curvature-induced 
terms are additionally suppressed in the long-wavelength 
limit, terms related to S arc neglected in what follows. 
It is then convenient to switch to Fourier space. 



u(r, f) = 



(3L 



E 



dqy 
27r ' 



-i{uJr 



where q = {qx,Qy), i^n = 2im/j3 with (3 = l/ksT are 
Matsubara frequencies, L = 2'kR is the circumference 

of the tube, and qx = n/R with integer n denotes the 
discrete transverse momenta. The action (2.1) then reads 
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2PL 



E 



dqy 
2-K 



u'^(-a;„, -q)A{hJn, g)u(a;„, q), 



(2.3) 



where A is given by the matrix 



BqxQy 



ml 



Bqxqy 



Mcol + Aql+fiql -t{B - ii)^ 



with A = B -h /i. 

Our strategy is to integrate out phonons, which is valid 
unless the electron-phonon coupling is very strong and/or 
soft phonon modes are present, see Ref. [11]. These re- 
quirements are met in our case. In order to proceed, we 
then need the inverse of -A, whose general form is straight- 
forward to obtain but not required here. The only trans- 
verse momentum relevant for the effective interaction is 
ga; = 0, since only the corresponding electronic states 
have to be kept at low energy scales. Therefore we fo- 
cus on this special limit henceforth. Putting qy = q and 
& = 0, the matrix A reads 



A = 



MljI + 

Mujl + l^q^ -i{B-iJi)q/R 

i{B-fi)q/R Ma;2+A/i?2 



with eigenvalues {i = 1,2,3) 

X.,=M[cjl + eUq)]. 
The inverse matrix is 



A-' 





V 





MujIr^+A i(B-ij,)q 

AaAa-Z?.^ A2A3H 

X2^3R A2A3 



Here ef{q) = v^q'^ with vt in Eq. (2.2), and 
el^{q) = {ujl/2)[l + q^R^TE{q)] 

with 

E{q) = [1 + 2q^R^{l - SB/x/A^) + q^R'^]^/^. 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



In the long- wavelength limit, qR<^l, further simplifica- 
tions arc possible. Although the theory is valid also for 
larger qR, we shall make this inessential approximation 
whenever appropriate. To lowest order in qR, we obtain 
for the three eigenmodes 
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e^{q) = vj.q , e^iq) = Vgq 
where 

vl = iB^/AM, 



A/MR^ 



(2.8) 



(2.9) 



The eigenmodes corresponding to Ai,2,3 are the twist- 
ing, stretching and breathing phonon modes, respec- 
tively. The sound velocities for the twisting and stretch- 
ing modes are vt and vs, and the breathing mode fre- 
quency is ujb- Prom Eq. (2.2) we obtain the following 
estimates 



Vt = 1.23 X 10^ m/sec, 
vs = 1.99 X 10^ m/sec, 

huiB = eVA. 
R 



(2.10) 



Note that both vs and vt are much smaller than the 
Fermi velocity, Vp = 8x 10^ m/sec. The stretching mode 
corresponds to the longitudinal acoustic (LA) phonon, 
while the twisting and breathing ones correspond to 
transverse acoustic (TA) phonons. 

Next we address the coupling of the phonon modes 
to electrons. The relevant electron states for a metallic 
SWNT arc spinors ippaa {v) [20] , where p = ± denotes the 
sublattice (the honeycomb lattice has two basis atoms), 
a = ± labels the two distinct K points, and a = ± the 
spin. Suppressing the a, a indices, the electron-phonon 
coupling corresponds to Hamiltonian densities H^^^ for 
the two K points. 
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Vi V2 



V2* Vi 



n 



K' 



*'e—ph 

Here the deformation potential is 



(2.11) 



(2.12) 



where gi « 16 eV in bulk graphite [28]. For a 2D graphite 
sheet this value should be multiplied by 3/2 [12], and 
we therefore take the estimate 511 rs 20 — 30 eV. The 
off-diagonal terms in sublattice space arise from a bond- 
length change, 



V2 = 526^'" (Uxx - Uyy + 2iUxy), 



(2.13) 



where 77 is the chiral angle, with 77 = for zig-zag tubes 
and ?7 = 7r/6 for armchair tubes. The respective coupling 
constant is (?2 ~ 1-5 eV [12]. Although both the free 
electron and the free phonon low-energy theory do not 
depend on chirality, the V2 term in the electron-phonon 
coupling does. 



III. INTERACTIONS FROM PHONON 
EXCHANGE 

We can now integrate out the phonons and obtain 
an effective action for the electrons that includes the 
phonon-mediated retarded interaction. To that purpose, 
we evaluate the Gaussian path integral for the phonons. 



j Puexpj-^ j [dq] 
+2V'^{-Un, -q)u{un, q)] } = e-^'^-^\ 
with the notation 



u^(-a;„, -q)A{uj„, g)u(a;„, q) 



(3.1) 



dq 
2^' 



(3.2) 



The result Sg-ph defines the phonon-mediated contribu- 
tion to the electronic action. In Eq. (3.1), we use the 
vector 



V(w„, q) = gi 



^Pa{<^n,q) 



M{q) 



(3.3) 



+ y 6-3*" [ Ji_ , q) + J2+ ,q)]N_iq). 

Fourier transforming Vi and V2 in Eqs. (2.12) and (2.13) 
specifies the auxiliary vectors M^(g) = {0,—iq,l/R), 
N^(g) = iq,iq,l/R), and N_iq) = N;(-g). Using 
fermionic Matsubara frequencies pn = (2n -|- l)7r/hf3 and 
analogous integration measure / [dp] , electronic densities 
are 



J[dp] 

Similarly, electronic currents are defined as 

'fpa{Pn,P)T±tpa{Pn + OJn,P + q)- 

(3.4) 

The fermionic (Grassmann) fields represent spinors, 
where we have suppressed the sublattice and spin in- 
dices, and Pauli matrices t± = ti ± zt2 act in sublattice 
space. The integration over phonons is now straightfor- 
ward, leading to the phonon-mediated action contribu- 
tion 

Se-ph = j [dq] V'^(-a;„,-g)A"^(a;„,g)V(a;„,q'). 

(3.5) 

This form will be analyzed in the remainder of the paper. 

To make further progress, it is convenient to invoke the 
Abelian bosonization scheme adapted to SWNTs [20]. 
To keep the paper self-contained, we briefly review the 
ingredients needed for our subsequent discussion. This 
requires first a unitary transformation in sublattice space, 
which transforms the free-electron Hamiltonian into the 
standard form containing right- and left-moving fields, 

V'pao- — ^r=R/L=± Upr'4'raai where 



u = 



_g-i7r/4 



%/2 



1 1 

i —i 



Under the unitary transformation, we have to replace 
T2 — > T3, Ta — > Ti, and ri T2. We then work in the 
transformed picture, and omit the tilde in the field oper- 
ators henceforth. The bosonization formula for the elec- 
tron operator is [20] 



Aaa{y) = 



: exp{i{rqFy + akpy + rV A-K(j)raa)} 



(3.6) 



where riraa are Klein factors, kp is the y-projection of 
the K point, and qp — \Ef\Ivf depends on the dop- 
ing level Ep. The chiral fields (f)raa are related to a 
pair of dual nonchiral boson fields, (paa = <t>Ra<T + 4>LoLa 
and 6acr = 4>Raa — 4>Laa- Switching to total/relative 
charge/spin fields, 4>c± = 5l]^(0ia ± 4>2a) and (j)s± = 
^ J2a- c^(0i<T i 02(7), and similarly for the dual fields, the 
electronic action including repulsive forward-scattering 
Coulomb interactions but without phonon exchange is 
described by a four-channel (A = c+,c—,s+,s—) Lut- 
tinger liquid [20], 



So = Y,U\dq] 



|^A(^n,g)|" 



{u;l + {uW). (3.7) 
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We ignore the tiny electron-electron backscattering pro- 
cesses [20] in this paper. The velocities are m° = vp/K^, 
where = 1 in all neutral channels (A ^ c+). For 
the charged channel, u'^,_^_ = vp/K^^, where the Lut- 
tinger liquid parameter K^_^ depends on the electron- 
electron interaction strength; K^_^_ = 1 for the nonin- 
teracting problem. In the presence of unscreened repul- 
sive interactions, K^_^_ ~ 0.2. In the rest of the paper 
we shall use the symbols u\,K\ without superscript to 
denote the values of the effective parameters renormal- 
ized by the phonon exchange. After this short review 
of bosonization for SWNTs, we now analyze the various 
terms in the phonon-mediated action S^-ph in Eq. (3.5) 
using bosonization. 



IV. ATTRACTIVE INTERACTIONS AND 
WENTZEL-BARDEEN SINGULARITY 

To begin with, there is a dominant forward-scattering 
term coupling to the total electronic density, which in 
bosonized form reads p{y) = {2/ y/n)dy(j}c+- The respec- 
tive contribution to Eq. (3.5) is 

'5/ = -^ [dq] p{-u)n,-q)Df{Wn,q)p{Un,q), (4.1) 
with the inverse propagator 

£>/(a;„,g) = M'^{-ujrr,-q)A-\uj,„q)M{iOn,q) 
Mwl{l + q^R^)+Aiiq^ 



This is the only phonon exchange contribution of order 
g\ . Due to the large deformation potential [as compared 
to the coupling 52], it is indeed expected to be the dom- 
inant one. 

Notice that the twisting mode, corresponding to the 
eigenvalue Ai, does not contribute to £>/. In order 
to make explicit the contributions from the other two 
phonon branches, we can rewrite D/ as follows, 



Df{un,q) 



1 



Z2((i,//.) 



Z-3{q,p) 



ME? 

where we use the functions 

i?2 



(4.2) 



Z2,3(«,A) =± 



AE{q) 



[-(l + g2E2)M43(g)+4Arj2]; 



E{q) is defined in Eq. (2.7). The contribution (4.1) thus 
only operates in the charged channel, there is no effect 
on neutral channels. Including the action in Eq. (3.7), 
the electronic action in the charged channel then takes 
the form 



1 

2vF 



S[ct>o+] = — I [dqU,+ {u;n,q)\'' ^1 + («°+) 



2 2 

r 



(4.4) 



^2(9,^) , Zsiq,)!) 



+ 4(1) 



where 



7 = 2vFgl/'K^nMR^. 



(4.5) 



The term related to Z2 is due to the stretching mode, 
while the term related to comes from the breathing 
mode. As pointed out above, the twisting mode drops out 
in this forward-scattering phonon exchange term, despite 
its importance for transport coefficients like the high- 
temperature conductivity [14]. Notice that the retarded 
phonon exchange contribution comes with a minus sign, 
indicating attractive interactions. 

Further simplifications in Eq. (4.4) are then possible 
in the long-wavelength limit. To lowest order in qR, 
Eq. (4.3) gives ^3(9, A) = 1 and .^2(9, A) = At;|gV(A - 
A)a;^, so that 



^2(9, m) _tJ'vlq^/Buj 



Zsiq,!^) 



Lot, 



vlq^ 



(4.6) 
(4.7) 



The forward scattering due to the stretching mode does 
not create a significant attractive interaction since it is 
suppressed by a factor [vs/vfY due to the linear dis- 
persion of this phonon mode [11]. A significant attrac- 
tive interaction can only be expected from the breathing 
mode. Since ujb corresponds to room temperature for 
SWNT radius i? « 4 A, see Eq. (2.10), we neglect retar- 
dation effects from the breathing mode in what follows. 
This approximation can be justified using the two-cutoff 
scheme of Ref. [19]. Putting uj'^+uj% ^ u)% in Eq. (4.7), 
interactions mediated by the breathing mode then sim- 
ply renormalize parameters in the effective LL picture. 
Concerning the stretching mode, see Eq. (4.6), we make 
no such approximation, and keep its influence on the WB 
singularity exactly. 

Straightforward algebra then leads to exactly the ac- 
tion studied before in Ref. [21], 



S[<p,+] = / [dq] 



I | 0c+(^«,g)|' 

2uc+Kc+ 



2 2 



h-'q' 



(4.3) with the parameter identiflcations 



i^+ = (t.°+)^-7M, 



(4.6 



(4.9) 



The corresponding LL parameter Kc+ is then given by 
Kc+ = Vf/uc+, leading directly to Eq. (1.1) with the 
radius Rb deflned as 
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2.9? 



T^hvpA 



2.4 ±0.9 A. 



(4.10) 



The uncertainty in the estimate for Rb largely results 
from uncertainties about the value for gi, see above. 
For this action (plus the standard free boson action for 
the neutral modes), all correlation functions have been 
worked out in Ref. [21], and one can simply adapt their 
results to the four-channel case encountered here. 

At first sight, when the radius R reaches the criti- 
cal value Rc = {K^_^)^Rb, a singularity is found, which 
corresponds to the WB singularity associated with the 
breathing phonon mode. However, due to the additional 
retarded interaction mediated by the stretching mode, 
i.e. the term ~ b"^ in Eq. (4.8), this singularity is reached 
at an even larger radius. Following Ref. [21], the condi- 
tion for the WB singularity is b/vs = Uc+, leading to a 
critical radius Rc = {K^_^)^Ro with 



Rq = {A/B)Rb = 3.6 ± 1.4 A. 



(4.11) 



The unexpectedly large value predicted for Rq suggests 

that it may be possible to reach this elusive singularity 
in practice. This may come as a surprise, since stan- 
dard reasoning tells us that a small ratio between sound 
and Fermi velocity implies negligible effects due to acous- 
tic phonon exchange, and this ratio is rather small in 
SWNTs. However, the effective phonon exchange at work 
here involves two phonon branches, namely the stretch- 
ing mode with £2(9) = vsIqI and the breathing mode, 
£3(5) = lob- These two modes reinforce each other in 
driving the system towards the WB singularity. In addi- 
tion, in SWNTs one has a very strong deformation poten- 
tial that provides very efficient electron-phonon coupling. 

An important conclusion in the context of nanotube 
superconductivity is that this forward-scattering phonon 
exchange process is very efficient in inducing attractive 
electron-electron interactions. To .judge the importance 
of the stretching mode, wo compare the scaling dimen- 
sions of various order parameters computed (i) under a 
non-retarded Luttingcr liquid description with LL pa- 
rameters w,c+ and Kc+ given in Eqs. (4.9) and (1.1), i.e., 
neglecting the stretching mode altogether, and (ii) keep- 
ing both the breathing mode and the retarded stretching 
mode, i.e. the action in Eq. (4.8), using a straightfor- 
ward generalization of exact results from Ref. [21]. Fo- 
cussing on the scaling dimensions for charge density wave 
(CDW) and singlet superconductivity (SC) order param- 
eters [20], this comparison is shown in Fig. 1. Remark- 
ably, concerning SC there is practically no difference at 
all. For the CDW order parameter, a significant differ- 
ence is observed only in the close vicinity of the WB sin- 
gularity. Sufficiently far away from the WB singularity, 
it is therefore possible to neglect the stretching mode and 
describe the electronic system as a Luttinger liquid with 
(possibly) attractive interactions. The stretching mode 
then only affects the location of the WB singularity but 



not the magnitude of the non-retarded attractive interac- 
tions. The LL interaction strength (including Coulomb 
interactions and phonon exchange) is determined by the 
parameter Kc+, which describes attractive interactions 
for Kc+ > 1. With Eq. (1.1), we provide an estimate for 
Kc+ as a function of the tube radius R, the interaction 
strength parameter K^^, and various material parame- 
ters. 



V. OTHER PHONON EXCHANGE 
CONTRIBUTIONS AND PEIERLS TRANSITION 

Let us now consider the remaining terms in the 
phonon-mediated action (3.5). First we combine the Ja± 
defined in Eq. (3.4), with a = 1, 2 = ± labehng the two 
K points, as Ji+ -I- J2- = iJc- + Jcdw and Ji_ -I- J2+ = 
-iJc- + JcDW- Here Jc- = S[dp]Y,a'^'^a'^'^'^<^ ^'^'^ 
Jcdw = /[(^J^] V'a'''iV'a. The reason for these defi- 
nitions becomes clear once we transform these operators 
back to real space and transform from sublattice to chiral 
space. In bosonized form, Jc-{y) = {^/V^) dyOc- and 

4 / 

JcDwiv) = — cos{2qFy + Vtt(Pc+) cos{^/t:4>c-) (5.1) 
Tra \ 

X s\n{\fK(j)s+) sin(-\/7r(/>s-) + [sin <r^ cos]^. 

These expressions show that Jc- is just the current 
density in the c— channel, whereas Jcdw is the inter- 
sublattice CDW order parameter [20]. Phonon interac- 
tions involving Jcdw are therefore important in relation 
to the Peierls transition and describe phonon backscat- 
tering processes. 

From Eq. (3.5), we then obtain several other contribu- 
tions in addition to the gi term in Sec. IV, namely S' 
and S" due to the (;| and gig2 terms, respectively. Let 
us start with S' = S'c_ + S'^£)viaj where 



S'c- = ~ J[dq]\Jc-{Wn,q)\' 

X [sin2(377)D+(a;„,g) +cos2(377)Z)_(g)] , 

and 

S'cDW = ~ / [^5] \ JcDwil^n,q)f 

X [cos'^{3r])D+{tOn,q)+sm^{3r])D_{q)] 
The inverse propagators are D-{q) = q^ /\\ and 

Mujl{l^q^R^)^ABq^ 



(5.2) 



MR^ 



q^R^) + 

\2\3R^ 

Z2{q,B) 



(5.3) 



(5.4) 



+ 



Zsiq^B) 



liq) 



The functions .^2,3 were defined in Eq. (4.3). In princi- 
ple, there is another contribution to S' from mixed terms 



6 



involving the product of Jc_ and Jcdw- However, such 
terms vanish for a doped SWNT by virtue of momen- 
tum conservation. In real space, this is reflected by the 
oscillatory cos{2qpy) factor appearing in Jcdw- Since 
this factor is absent in Jc_, we ignore the mixed terms 
in what follows. The second contribution to the action 
reads 



S" 



■ sin(3r?) j [dq] 



(5.5) 



X pc+{-u}n,-q)D^{un,q)Jc-{oJn,q), 



with 



Dm{uJn,q) = 



(5.6) 



1 



Z2{q) 



+ 



^l + 4{q) u^l + eliq) 



where ^2,3(9) = t(1 - q^R'^)MR^el.^{q)/AE{q), with 
E{q) given in Eq. (2.7). At long wavelengths, this yields 
Z2{q) = -vlq^/ujl and Zz{q) = 1. Then Eq. (5.6) sim- 
plifies to 



Dm{tOn,q^R~^) = 



(5.7) 



Again, there is in principle another contribution that 
arises from mixing pc+ with Jcdw- However, this contri- 
bution also vanishes by virtue of momentum conservation 
and is disregarded. 

Next, we analyze these additional terms using 
bosonization, starting with the forward-scattering term 
S'^_, which can be treated along the same lines as the 
forward-scattering term in the c+ channel, see Sec. IV. 
Neglecting retardation in the breathing mode as above, 
the long-wavelength action for the bosonic field 9c- in- 
cluding is 



S[Gc-] = ^ J[dq]\0c-{uJn,q)\^{i^. 



sin^(3?7) + 



L0% 



UJt, 



+ cos2(37?) 



q^R^ 



2 „2 



with the parameter 7' = 2wFg|/7r2?j,Mi?^. Amusingly, 
the contributions of different phonon modes are weighted 
by chirality-dependent factors. In particular, for zig- 
zag tubes (?7 = 0) only the twisting mode contributes, 
whereas for armchair tubes (77 = 7r/6) only breathing 
and stretching modes do. Again we obtain a WB sin- 
gularity but now in the c— channel at a corresponding 
critical radius R'q = 2g2/7T'^UvFp. Inserting the parame- 
ter values given above, due to the small ratio 52/51, the 
critical radius is exceedingly small, R'q « 0.02 A. Using 



the same reasoning as in Sec. IV, the effect of S'^_ again 
leads to a LL action in the c— channel, with a renormal- 
ization of Wc- and Kc- - However, such renormalizations 
amount to tiny changes of the order 10^"^ even for very 
thin SWNTs. Therefore S'^_ does not imply observable 
consequences and is omitted henceforth. 

Next we focus on S" , which mixes the c+ and c— sec- 
tors and constitutes another forward scattering mecha- 
nism. First, we observe that this term vanishes for zig- 
zag tubes (77 = 0). Neglecting again retardation in the 
breathing mode and keeping only terms to lowest order 
in qR, we find from Eq. (5.7) 

S" = eK [[dq]cj)c+{-uJn,-q)0c-{oJn,q) ^"^^2 2 ' (5-9) 

with e = (72/51 and k, = sm{3ri)vFRB /fiR- Introducing 

the notation <f>^ = {(j)c+,0c+,(l)c-,9c-)- the total action 
[including Eq. (4.8) and S"] can be written as 

S=ll[dq]^'^[S^ + eS,]^, (5.10) 

where is a block-diagonal matrix in (c-f,c— ) space, 
and S_i contains the mixing term due to S". Since e <C 1, 
we proceed by expanding to first order in e. This allows 
the explicit calculation of scaling dimensions of all or- 
der parameters of interest. Details about this calculation 
are contained in an appendix, and the main results are 
shown in Fig. 2. Clearly, the effect of S" is also very tiny, 
even for very thin tubes. Neglecting this term, only the 
term S'^j^^^r describing backscattering phonon exchange 
remains. 

Let us then study the S'^^^y term. As this describes an 
effective interaction involving Jcdw, this part monitors 
the possibility of a Peierls transition. In order to keep 
the discussion transparent, we focus on the static limit. 
Then Eq. (5.3) reduces to 

S'cDW = j [dq] \ JcDw{(^n,q)f- (5.11) 

Transforming to real space, using the bosonized expres- 
sion of Jcdw in EQ- (5.1), we find a term involving Aqp 
oscillatory contributions that violate momentum con- 
servation and average to zero, and in addition a tiny 
renormalization of the LL parameters u\ and Kx. This 
renormalization represents a direct generalization of the 
one/two-channel theory of Voit and Schulz [11] to the 
four-channel case of interest here. We do not give any 
details here, since the relevant dimensionless coupling 
constant is Rq/tt'^R, and the resulting renormalizations 
are again of order 10~^, i.e. unobservable. Neverthe- 
less, we confirm the finding of Ref. [15] that the relevant 
Peierls distortion mode for armchair SWNTs is given by 
a twist distortion, see Eq. (5.3). For arbitrary chiral an- 
gle 7], however, the relevant distortion involves a linear 
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combination of all three phonon modes. Finally, we men- 
tion that on half-filling, the oscillatory contribution could 
survive. However, its scaling dimension suggests that 
this term is strongly irrelevant unless electron-electron 
backscattering gaps out the neutral channels. For very 
thin tubes, these gaps may be sufficiently pronounced [20] 
to allow for Eq. (5.11) to be a relevant perturbation. Pro- 
vided ifc+ < 2, one could possibly expect a small (twist) 
Peierls gap in this nearly commensurate situation, see 
Rcf. [15]. However, this would require a detailed analysis 
of additional phonon exchange processes neglected above 
due to incommensurability. Away from half-filling, we do 
not see a possibility for this gap to survive. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied acoustic phonon ex- 
change mechanisms operating in an individual metallic 
SWNT. The description of phonons within an elastic con- 
tinuum model allows for explicit analytical progress, and 
we have shown that integrating out the phonon system 
can produce very substantial attractive electron-electron 
interactions. For sufficiently large tube radius, the sys- 
tem can be described by a standard Luttinger liquid the- 
ory, where the interaction parameter Kc+ is renormalized 
away from the value K^j^ < 1 describing the theory with- 
out phonons but fully taking into account Coulomb inter- 
actions. We have shown that Kc+ > 1 corresponding to 
attractive interactions is possible, in particular once the 
Coulomb interactions are externally screened. These at- 
tractive interactions are basically mediated via a breath- 
ing phonon mode. For thin SWNTs, one also has to 
take into account a retarded branch describing a stretch- 
ing mode. The combined effect of the stretching and 
the breathing mode then produces a Wentzel-Bardeen 
singularity at a critical radius Rq « 3.6 A. Nanotubes 
may then offer the possibility to reach this singularity 
in experiments. Remarkably, all other mechanisms, e.g. 
backscattering phonon exchange, turn out to be remark- 
able inefficient in mediating effective electron-electron in- 
teractions. In particular, based on our analysis, it should 
be extraordinarily difficult to observe traces of the Peierls 
transition. As our treatment is essentially exact, we be- 
lieve that this also invalidates some earlier theories that 
were mainly based on mean-field approximations. 

As mentioned in the introduction, the phonon- 
mcdiatcd attractive interaction can only drive the sys- 
tem into a true superconducting phase if inter-tube cou- 
plings are taken into account. This will be studied in de- 
tail elsewhere [23]. In order to gain preliminary insights 
about inter-tube phonon exchange effects in a rope of 
SWNTs or a multi-wall nanotube, let us briefly contem- 
plate a simplified model of two parallel SWNTs and study 
phonon exchange in the presence of elastic couplings be- 
tween both tubes. We shall focus on the breathing mode 



alone, which will produce the most important effects. In 
the long- wavelength limit, it can be approximated as a 
dispersionless mode, leading to the free-phonon action 

S=^f[dq\ ^ Ui{-u;n,-q)Aij{uJn)uj{uJn,q), 

where Ui is the radial component (uz) of the displacement 
field on tube i. 



-atuj% ojl + oj%{l + ai) 



and at,ai are dimensionless parameters. While one 
can compute these parameters from a microscopic force- 
constant model, here we take them simply as phenomeno- 
logical parameters. On general grounds, they must obey 

< at, I < 1, where at i = in the absence of inter- 
tube couplings. With total electronic densities pi on tube 

1 = 1,2, the dominant electron- phonon interaction is due 
to the deformation potential (3.3), producing a contribu- 
tion 

{gx/iirB?) I [dq] ^ pi{-Wn,-q)ui{u;n,q) 

to the action. We can diagonalize this action by in- 
troducing symmetric and antisymmetric combinations, 
Us/a = {ui ± U2)/\/2 and ps/a = {pi ± P2)/V2. Integrat- 
ing out the phonons, the forward-scattering contribution 
in the c+ channel leads to a Luttinger liquid, but with 
renormalized parameters. 



c+\ 



j=s,a 



K + {uj,c+fq^^ 



(6.1) 



where Us/a,c+ = ^f'/^s/a,c-i-- The interaction parameters 
are 



K 



(6.2) 



where f^/a = 1 + aiT at- 

There are now two types of superconducting order pa- 
rameters. Focusing on the most important singlet super- 
conducting fluctuations, there is one operator describing 
a Cooper pair on the same tube, O 



SSCO' 



and another 

one where the electrons forming the Cooper pair reside 
on different tubes, Ogg^Q, with explicit form 

^SSCO ~ Cri>l,raai>l,-r-a-CT, 
raa 

^SSCO ~ ^ CrV'l,ra<TV'2,-r-a-<T, 
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where "01/2, raa is the electron spinor on tube 1/2. From 
Eq. (6.1), the scahng dimensions of these operators fohow 
as 



^sTco = - I 6 + 



A inter 

^ssco 



1 



1 



6 



K 



K, 



a,c+ 



(6.3) 
(6.4) 



Two conclusions can be drawn from the inspection of 
these scaling dimensions. First, for inter-tube couplings 
with at sufficiently larger than a; , attractive interactions 
within a given tube are increased. Second, whenever 
Ka,c+ > 1, which holds true in the absence of Coulomb 
interactions, the intra-tube superconducting order pa- 
rameter is always more relevant than the inter-tube one. 
In other words, we expect that Cooper pairs predom- 
inantly form on a given tube but not across different 
tubes. Remarkably, this conclusion does not depend on 
the values of a;,*. Clearly, this finding will be impor- 
tant in the theoretical analysis of experiments on SWNT 
ropes and related systems. 

Finally, we briefly address the role of electron-electron 
backscattering that has been neglected in our study. 
From Ref. [20], backscattering is expected to open tiny 
gaps that depend exponentially on the radius R. While 
for i? > 5 A, the gap is extremely small and of no rele- 
vance in practice, for ultrathin tubes (i? « 2 A) the gaps 
may be more important. However, in this limit also the 
bandstructurc may change in a profound way duo to hy- 
bridization of TT and a orbitals, and the analysis of Ref. 
[20] cannot simply be taken over. In addition, the pres- 
ence of the Wcntzol-Bardccn singiilarity poses an intrin- 
sic lower limit to the validity and self-consistency of our 
eff'ective low-energy approach. With the predicted value 
Ro w 3.6 A, wc therefore do not expect that backscatter- 
ing leads to dramatic changes to the picture put forth in 
this paper. 
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APPENDIX A: PERTURBATIVE TREATMENT 
OF S' 

In this appendix, we briefly describe some details con- 
cerning the perturbative treatment of S", see Eq. (5.9), 
in Sec. V. With the notation = {(l)c+,0c+,4'c-,()c-), 
the full action in the total/relative charge channels 
is Eq. (5.10) with the block-diagonal matrix 5q = 
diag(A^+, A^_), where 



"c+ _2 _ q /VF 




VfQ 



vpq 



-ILOnq 



Furthermore, the mixing term leads to a matrix S_i, 
whose only non-zero entry is 

(S;i)i4 ^ (S_-i)ii = Ku:lq^ l{ujl+vlq^), 

where k, is defined immediately after Eq. (5.9). We then 
expand to first order in e. 



(^0 + e5i)-i = S^' - eSo'S.So' - 

H 



q-l 







Some algebra yields for the matrix H_ the following result: 

H- 

where 



{iVFUJn/q)'P vjpV 



v 



V 



E 



/3 



v 



{iVFi^nlqjV 



E 

0=±, 



2 „2- 



Here i = 0, 1, 2 and we denoted vi± = vs, V2± = vf- The 
velocities vo± = v± and the coeflacients Ff}, Cp are given 
in Ref. [21]: 



2v\ 



Fa = 



Uc+ Vj - ' 



C0 . . 

Finally, the coeflacients Qf are defined by 

Qt 



-4 


(4 










-4)(4-4) 








-vl){vj,-vl) 



With this, it is straightforward to compute correlation 
functions for the bosonic field. The T = equal-time 
correlation functions ($j(a;)$j(0)) are given by 

($.(x)$,(0)) = -^ln([x2+a2]/L§), 

where a and Lq are UV and IR cutoff lengths. The only 
non-zero entries of Gij are G33 = G44 = 1 and 
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Gil = '^vfFi}/v0, 



Gi4 = -cKVp^^FfjQl /v^jj, 

G22 = ^Uc+Cij/VF, 



0,i 

With this, the scahng dimension of any operator 
exp[ia^^{x)] is found as a^Ga/A-K. 
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FIG. 1. Comparison of the scaling dimensions as a function 
of radius R (in A) for SC (lower curves) and CDW (upper 
curves) order parameters. Here we take the deformation po- 
tential coupling constant gi — 25 eV. The dashed curves give 
the non-retarded LL prediction, while the solid curves rely on 
Eq. (4.8). The dotted vertical line denotes the location of the 
WB singularity. In the case of the SC order parameter, the 
dashed and the solid curves cannot be distinguished. 
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1.4 
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FIG. 2. Comparison of the scaling dimensions as a function 
of radius R (inA) for the SC (lower curve) and the CDW 
(upper curve) order parameter, computed with e = (dashed 
curves) and e = 0.06 (solid curves). The calculation employs 
Eq. (4.8). Taking gi = 25 eV, the value e = 0.06 corresponds 
to 32 ~ 1-5 eV. We consider the case of an armchair tube 
{ri — 7r/6), where the effect of S" is most pronounced. 
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